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The concept of mesons modified by the nuclear medium has been used for many years 
]]]]. It has been generally assumed that although the properties of pions and kaons (mass, 
decay constant, coupling constants, etc.) are changed by the presence of the medium, their 
pseudo-Goldstone nature is preserved.[] In the vacuum, the masses of pions or kaons are 
known to be proportional to the square root of the current quark masses, according to the 
Gell-Mann-Oakes-Renner (GMOR) relation ||. However, in nuclear matter this need not 
be the case. 

In this paper we show, that the nature of charged pionic excitations in a dense isospin- 
asymmetric medium, or kaonic excitations in a non-strange medium, is radically different 
than what is seen in the vacuum results. With no dynamical assumptions other than the 
fact the current quark masses are sufficiently light, and that expectation values of physical 
operators do not diverge, one can show that in such a medium there must exist at least one 
pseudo-Goldstone mode with an energy that scales with the average current quark mass as 
m d where d > |. This is in sharp contrast with the case of the vacuum or isoscalar matter, 
where the energy of the Goldstone modes goes as m 1 / 2 . 

Moreover, if one makes highly plausible assumptions about the behavior of the chemical 
potential for the system as one varies the quark masses, then one can show even more peculiar 
behavior. In the vacuum one has pairs of charged pseudo-Goldstone excitations (7r + -7r~ 
or K + -K~). As one goes to a dense isospin-asymmetric medium (such as dense neutron 
matter), however, given these assumptions one can show that only one member of this pair 
survives as a pseudo-Goldstone mode. In addition, the excitation of the surviving pseudo- 
Goldstone boson is proportional to the current quark mass m q itself, rather than to ^m 9 , 
as in the vacuum case. Specifically, we find that in the chiral limit in neutron matter there 
exist 1) a pseudo-Goldstone mode with quantum numbers tt + with excitation energy ~ m q , 



A mode can be identified as a pseudo-Goldstone mode, if its excitation energy vanishes in the 
chiral limit of zero quark mass. 
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2) pseudo-Goldstone mode with quantum numbers tt° with excitation energy ~ ^ m g, and 

3) there are no pseudo-Goldstone modes with quantum number of tt~ . 

Our analysis is similar in many respects to that of Ref. 0|, however that work implicitly 
assumes that the medium has low density. Mathematically, this equivalent to expanding in 
density before expanding in the quark mass. This is clearly appropriate for mesonic atoms, 
where the meson feels only the tail of the nucleus at a fraction of the nuclear saturation 
density, po- Here we are interested in studying sufficiently dense systems so that it becomes 
appropriate to consider chiral limit with the density kept constant. We will present some 
model estimates indicating that this situation may be realized in nature if the density of 
matter is of the order of a few po- We will also demonstrate explicitly for a simple model 
that the chiral and low (isovector) density limits do not commute: an expansion in the small 
parameter associated with one limit is singular in the expansion parameter associated with 
the other limit. Thus, physics associated with the two limits will be quite different. 

The paper is organized as follows: First, we derive our general result of existence of 
unusually soft pseudo-Goldstone modes in nonsymmetric nuclear medium. Next, we present 
a simple model illustrating that behavior and give some numerical estimates of where we 
might see these unusually soft modes. Finally, we discuss the case of kaons in non-strange 
matter. 

We find it useful to introduce the function dim x (X) = linin^o (m ^pb- ) , which we call 
the chiral dimension of the quantity X. If near the chiral limit X scales with the quark 
mass as a power law in m, i.e. if X ~ m a , then function dim x (X) extracts the power a. 
For instance, in the vacuum dim x (m„-) = 1/2 and dim x (F 7r ) = 0. The result is not modified 
by the presence of logarithms i.e. dim x (m a log ra m) = a. 

The starting point of our analysis is similar as in the derivation of the standard GMOR 
relation. We examine the expression for the double commutator of the QCD hamiltonian 
density at point x = 0, 7i(0), with the axial vector charges Q% — j d 3 xJ£ (x), where 
is the axial vector current. This expression can be evaluated directly through the use of 
commutation rules for the quark fields: 
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[Qt, {Ql H(0)]] = m {X a /2, {A72, M}} ^(0) , any a = 1, ... 
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(1) 



where A a are the Gell-Mann flavor matrices, M = diag(m M , rrid, m s ) is the quark mass matrix, 
and curly brackets denote anti-commutators. Next, we take the matrix element of both sides 
of the above operator relation in some physical state |C). In all of our applications this state 



linear combination of quark condensates of various flavors. For simplicity of notation, let 
us consider the case of two flavors, with m u = rrid = m. The generalization to 3 flavors is 
straightforward. We then obtain 



(C\ [Ql [Ql H(0)}} \C) = m(C\ (uu(0) + dd(0)) \C) = m{qq)c , a = 1, 2, 3. (2) 



On the other hand, the LHS of the preceding equation may be evaluated by inserting a 
complete set of intermediate states, yielding a sum rule. 

Since the system \C) is spatially uniform, these intermediate states can be labeled by their 
three-momentum, and some additional label j. We denote them as \j,p)- It is convenient to 
measure the three- momentum, p, and the energy, Ej, of the intermediate states relative to 
the ground state \C) at rest. The total energy and the total three- momentum of the state 
\j, p) and of the ground state \C) both form Lorentz four- vectors. Their difference is therefore 
also a four-vector, and the unity can be decomposed in the following Lorentz-invariant way: 



Inserting Eq. (|3|) inside the LHS of Eq. ([[J), and using the fact that the ground state 
and the intermediate states are eigenstates of the hamiltonian H = J d 3 xTi,(x), namely 
H\C) — Ec\C) and H\j,p) = (Eq + Ej)\j,p), we obtain the following relation: 



will be taken to be spatially uniform. The matrix element of the RHS of Eq. ([D) yields a 




(3) 



E 



(c\Jto(oMp 




m(qq)c , 



(4) 
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We may use the following operator relation 



D a (0) = V(0)i7 5 r a V(0) 



(5) 
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to relate the matrix elements of J5 to the matrix elements of D a . We find 

777 

(C\ J 5 y0)|j,p = 0) = -i —(C\D a (0)\j, p = 0). (6) 

^3 

Using Eq. (|^) we may rewrite Eq. in the equivalent form which is familiar from the 
derivation of the GMOR relation: 

E H^r \(C\D a (0)\j,p = 0)| 2 = -(qq) c . (7) 
j i o\ J 

Let us now recall the usual arguments leading from the preceding general result to the 
GMOR relation. In the GMOR case, the state |C) is simply the vacuum. First, we note that 
all excitation energies Ej are positive, since the vacuum is, be definition, the lowest-energy 
state of the system. Therefore all components of the LHS are positive. Next, consider taking 
the limit m — > on both sides of Eq. (^). The RHS is assumed to go to a nonzero constant as 
a result of the spontaneous chiral symmetry breaking, i.e. dim x ((gg)) = 0. On the LHS only 
the terms with pseudo-Goldstones can contribute. This is because dim x ((0\D a (0)\j,p = 0) = 

0, and thus only the states for which dim x (Ej) = 1/2 can contribute to the sum rule in the 
chiral limit. In the vacuum, Ej is just the mass of the excited particle (pion). Furthermore, 
from isospin symmetry we have m n o = m v + = m v ~ . 

In the present case our state |C) is not the vacuum. Rather, we take |C) to be the 
ground state of the system subject to a constraint which requires that the expectation value 
of some local (space-independent) operator c(x) (which commutes with H) has a fixed value, 

1. e. (C\c{x)\C) = c = const. In our case we are interested in medium with a given flavor 
density, and c(x) = fi u ii + u + Hdd + d or c(x) = 3fiB(u + u + d + d) + jii = \{u + u — d + d), where 
/x's are the corresponding chemical potentials for up and down flavors, or alternatively for 
the baryon number density and the isovector density operator u + u — d + d. If c ^ 0, then we 
have no guarantee that Ej is positive definite. The operator Jg may connect to states with 
a lower energy. 

There are two distinct cases to consider: 1) the operator c(x) commutes with the axial 
vector charges or 2) the operator c(x) does not commute with the axial vector charges. As 
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already noted in Ref. [[|, if c(x) commutes with Q§ (or equivalently with D a ), then the states 
with negative Ej do not contribute in the sum rule. This is simply because in this case D a 
connects |C) only to intermediate states with the same value of c, and, by definition, \C) is 
the lowest energy state of a given value of c. 

One physical situation where the constraint commutes with the D a operators is the case 
of finite baryon density with zero isovector density (isospin symmetric nuclear matter). In 
this case the conclusions from the sum rule (0) are obtained analogously to the vacuum case, 
and, so long as (qq)c 7^ in the chiral limit, we conclude there must exist pseudo-Goldstone 
excitations with Ej ~ with quantum numbers of 7r°, tt + and 7r~.Q Isospin symmetry 
causes the neutral and charged excitations to be degenerate in energy. Thus, from the point 
of view of the GMOR relation, symmetric nuclear matter behaves similarly to the vacuum. 

Now we come to the main topic of our considerations. Take nonsymmetric nuclear 
medium, i.e. pi=\ = (C\{u + u — d + d)\C) 7^ 0. The neutral pseudo-Goldstone excitation still 
behaves in the usual way, since the neutral axial vector charge commutes with the third 
isospin component of the isospin charge, Q\. Hence, as already remarked in even in 
nonsymmetric matter we have a neutral pseudo-Goldstone excitation such that E^o ~ \/m. 
The case of charged pionic excitations, however, is radically different for two reasons. The 
first is that the isovector constraint does not commute with the charged axial vector charges 
Q\ and Q\. The second difference is related to the existence of a second sum rule which is 
trivially zero for an isoscalar medium but not for a medium with nonzero isospin density. 
This sum rule is derived rather easily. Using the fact that Q§, J^ (0) = ze afcc Jo(0), where 
J° is the vector current, and inserting Eq. (|3|) inside the commutator, we obtain 

2 2 

ffi Tfi 

Pi=i = T,^r^\{3-^D~mC)\ 2 -Y.T^r^\{3 + MD + mC)\ 2 , (8) 
where we have decomposed the sum over j into two classes of states: those with isospin one 



2 The case where (qq)c = is interesting in its own right. This problem is discussed in detain in 
ref. 1 
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unit more or less than in \C). Equation (0) may be decomposed similarly, giving 

-(gg>c = E 2 .R, |r, \^P\D-mC)\ 2 + Y, 2 \ E . lE . \(j + ,p\D + mC)\ 2 . (9) 

Now consider the sum rule given in Eq. (§) as one approaches the chiral limit. In isovector 
matter, the left hand side of this sum rule is nonzero, by definition. Accordingly, at least 
one term in this sum rule must be nonzero. However, all terms in the sum rule go as mj Ej 
times matrix elements. Since the matrix elements are assumed to be finite in the chiral 
limit, one sees that the only way any term can make a contribution in the chiral limit is if 
dim x (.Ej) > 2/3. This completes the demonstration of the first point of this Letter. Note, 
that this relation is written as inequality rather than an equality since it is possible that the 
chiral dimension of the matrix element could well be > 0. Indeed, as we will show below, 
given reasonable assumptions it is in fact > 0. 

To proceed further we need to make the assumption that the chiral dimension of the 
isovector chemical potential dim x (/xj = i) = 0. Since we are implicitly taking the chiral limit 
at a fixed isospin density pi=i, we can treat pi=\ as an external parameter ; independent 
of the chiral parameter — the assumption which we must make is that the difference in 
energy density between this state and the the lowest isospin symmetric state is independent 
of m. Although we cannot prove that dim x (/// = i) = from first principles, we can present 
the following physical argument in its favor: The isovector interaction has a contribution 
produced by the p-meson exchange. Suppose we place an object of isospin J 3 in the isovector 
medium. The interaction is {gp/m 2 ) ^7=1/3, and the corresponding chemical potential is 
{g 2 p /rn 2 p ) pi=\. Here g p and m p are p-meson coupling constant and mass in the medium. 
Their chiral dimension in the vacuum is 0, and, unless something very unusual happens, this 
is also true in medium. Therefore p-exchange produces dim x (p,/ =1 ) = 0. It is very unlikely 
that this result could be altered by other processes — they would have to exactly cancel the 
p-exchange mechanism. 

Equipped with the assumption that dim x (p/ = i) = 0, we can continue the analysis of the 
sum rules (§) and @. Suppose for definiteness that pi=\ > (as it is for neutron matter). 
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From definition of the chemical potential as the minimum energy needed to lower the isospin 
by one unit, we have 

E jT = fi I=1 ± 6E U , with 5E U > . (10) 

Therefore the energies of states with one unit of isospin less than \C), Ej_, can never go to 
zero. In the chiral limit only the states with Ej — > can contribute to the sum rules @ 
and (P), and as a result only the states j + contribute. Note, this is very different from the 
vacuum case — it shows that in neutron matter the ir + modes can be pseudo-Goldstone 
modes while the 7r~ cannot. Restricting the sum rules to the j + modes only gives 

pi=x = ^{^T,^j\^P\D + mc)\^ , (ii) 

and 

hm {{qq) c ) = -hm ™ \(j + ,0\D + (O)\C)A . (12) 

m-0 m-*0 y j+ 2\E j+ \E J+ J 

Now, the sum rule fllT]) contains only semi-negative contributions. Comparing the chiral 
dimensions on both sides, it follows that there must exist a mode for which 

= 2-3 dim x (£ J+ ) + 2 dim x ((j+,P = 0\D + (0)\C}) . (13) 

The sum rule ( |T2"D may contain both positive and negative contributions, since the sign of 
Ej + is not restricted. That means that in principle there may be cancellations of the leading 
chiral dimension on the RHS of Eq. ([12]), and we obtain the following inequality: 

dim x (gg) > 1 - 2 dim x (E j+ ) + 2 dim x ((j+,p = 0|£>+(0)|C» . (14) 

The conjunction of conditions (p~3|) and (Hj) gives dim x (Ej + ) < 1 + dim x ((qq)) c ■ The 
inequality becomes equality unless there is an exactly cancellation of the leading order 
contribution on the RHS of Eq. ([12]). For instance, the equality is automatically the case if 
there is only one state \j+,p = 0) which becomes a pseudo-Goldstone mode. Also, potential 
cancellations are not associated with any symmetry, hence it is difficult to imagine that the 
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leading chiral dimension can indeed be exact canceled on the RHS of Eq. fll2|). With no 
cancellations we have 

&m x (E j+ ) = 1 + dim x ((qq) c ) . (15) 

If we are in the spontaneously broken phase, then we expect that dim((qq) c ) = 0, as in 
the vacuum. Clearly, this dimension has at least to be positive semi-definite, if the chiral 
condensate is to be well-behaved in the chiral limit. This leads to the following final result: 

dim x (£ J+ ) = 1 , dim x ( J D+) = 1/2 , dim x (J+) = 1/2 , (16) 

where the third equality follows from Eq. @. Note that this behavior is radically different 
than in the vacuum, where we have 

dim x (m w ) = 1/2 , dim x ( J D + ) = 0, dim x (J 5 +) = 1/2 . (17) 

Of the three pionic pseudo-Goldstone modes in the vacuum, only only 2 survive as pseudo- 
Goldstone modes in dense isovector matter. For negative pi = \ the negative charge excitation 
disappeared, the positive charge excitation became unusually soft, E n + ~ m, and the neutral 
excitation retained its chiral dimension, E^o ~ \fwi. If pj = \ > 0, then the positive and 
negative excitation change the roles. 

The fact that only one member of the pair of charged pseudo-Goldstone mesons remains a 
pseudo-Goldstone mode in dense isospin-asymmetric matter is quite natural. The medium 
itself breaks the symmetry between the members of the pair — thus one would be quite 
surprised if they had the same energy. One expects the symmetry breaking induced the 
medium to split the degenercy between the two states. However, as one approaches the 
chiral limit a pseudo-Goldstone mode must go to zero excitation energy. Thus, if both 
modes were pseudo-Goldstone modes they would both have to approach zero and hence 
would become degenerate. 

Before giving an illustrative example, let us recapitulate our derivation, listing all as- 
sumptions made on the way. First, let us stress that what we were after is a fundamen- 
tal result derived without any explicit reference to dynamics, microscopic structure of the 
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modes, etc. In matter these modes are undoubtedly quite complicated, involving particle- 
hole excitations, etc. We have not made any dynamical assumptions. These were our key 
ingredients: 

1. The chiral limit is taken first, while the isovector density is kept constant. 

2. "Reasonableness" conditions are that expectation values of operators, e.g. (qq)c, do 
not diverge in the chiral limit. If this were not true, than the chiral limit would not 
make sense in the isovector matter. 

3. The isovector chemical potential is assumed to be nonzero in the chiral limit. 

4. We also assumed that there are no exact cancellations of the leading chiral powers in 
the sum rule (^|). This assumptions is equivalent to having no (accidental) cancellations 
from a priori possible multiple branches of pseudo-Goldstone modes. If there is only 
one such branch, this result follows trivially. 

Assumptions 1) and 2) are sufficient to show that there exists a pseudo-Goldstone mode for 
which dim x (Ej) > 2/3. Assumptions 1) - 4) give dim x (Ej) = 1 for this mode. 

Now, let us present a simple model which will illustrate the behavior of pseudo-Goldstone 
modes in nonsymmetric medium. Consider the pion propagating in an isovector medium 
and interacting with it through the p-meson exchange. Moreover, in this toy model, we will 
assume that the p-meson exchange is the only interaction between the pions and the medium. 
As is well known ||, if all vector-isovector interactions are mediated by the p meson, then 
consistency of the p exchange picture with the various soft-pion theorems requires the KSFR 
H relation m 2 p / ' g 2 = 2F 2 to be exact. The inverse propagator for the mesons with four- 
momentum p is (G^p) 1 * 1 ) -1 = p 2 — m 2 =ppoA, where A = (g 2 /m 2 )pj = i. The (positive energy) 
poles of G{p) for the pion at rest (p = 0) occur at 



E + = -A/2 + ^A 2 /A + ml 




(18) 
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In the chiral limit, E + — > m\jA — this is the unusually soft pseudo-Goldstone mode, 
with dim(E + ) = 1, and E~ —>■ A — this is no longer a pseudo-Goldstone mode, since 
dim(E~) = 0. This is precisely the behavior which we predicted. Obviously, the neutral 
pion in unchanged by the isovector interactions, and dim x (E + ) = 1/2. 

Let us also examine the matrix elements of the operator D 1 * 1 . We have, from definition 
(H), D ± = d^J^/m = (F w m 2 ,/m) 7r ± , where ^ is the (charged) pion interpolating field, 
and F n is the pion decay constant in the vacuum. Let abbreviate the states corresponding 
to excitations E ± by |±). Using | < =1= | tt^ | C> | 2 = \im po ^ E j 2E j (p - E j )G ± (p ,p = 0), where 
j = + or — , we find that in the chiral limit |(+|£) + |C)| 2 — > 2F 2 m\j \A 2 m 2 ) ~ m, in 
compliance to the general result (|13|) , and | (— \D + \C) | 2 ~ 1. Only the positive charge mode 
saturates the sum rules (H) and @. Explicitly, we get 

PI=1 = F 2 A = 2( 9 -^) 2 p I=1 , (19) 

and 

F 2 2 

Km {qq)c = . (20) 

Equation (|TSp means that, as expected, consistency in the chiral limit requires the KSFR 
relation to be exact. Relation fl2"U| ) also shows the formal consistency, since exactly the same 
equation is obtained by considering the sum rule for the neutral pion. Since the neutral pion 
is not affected by the p-exchange, we immediately get Eq. (^0|). 

Under what circumstances is the analysis here useful? - - i.e. at what isospin density 
does the system go from effectively having three pseudo-Goldstone modes whose energies 
all have chiral dimension 1/2 to having two pseudo-Goldstone modes whose energies have 
chiral dimension of 1 and 1/2. Formally, the analysis is valid for any nonzero isospin density 
provided we go to the m q —>■ limit. On the other hand, in nuclear physics the value of m q , 
though small, is nonzero. Thus, the question of interest is under what circumstances is the 
quark mass small enough so that the results of m q — > limit are close to the physical results. 
The toy model considered above, although probably quite unrealistic, gives us considerable 
insight. 
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Consider Eq. (0); one approaches the limit in which the present analysis applies when 
\A\/2 ^> m n . We note that in this toy model A = (pi=ig 2 p ) /m 2 p is the chemical potential 
associated with isospin, pi=\. Also in this simple model, is the energy of the pseudo- 
Goldstone mode when pi=\ = 0. In the general case it may be different from m n , so let 
us denote it by E . Thus, a reasonable general criterion for the region of applicability 
of our analysis is when /ij=i 3> Eq. In our toy model, this condition is equivalent to 
|/0j=i | 3> p = 4:F 2 m n ~ 3.7po, where po is the nuclear saturation density. The above estimate 
used the vacuum values of F n and m n . However, this probably leads to an overestimation 
of p, because matter with large isovector density also has large baryon density, ps, which 
considerably reduces the value of F v from its vacuum value In turn, p should be 

substantially reduced by the presence of ps- 

We also note one more fact illustrated by our toy model. It is apparent that the chiral 
limit and the low-isovector-density limit do not commute. In our model in the chiral limit 
the expansion parameter is a = m^F 7r /p/ =1 - - this is singular in the isovector density. 
Conversely, in in the low density limit the expansion parameter is 1/a, and this, it turn, is 
singular in m n . We believe this is a manifestation of a general result. 

Now let is turn our attention to kaons, since recently the possibility of 5*-wave kaon 
condensation has been extensively discussed |l5|- p0| , p!4| . The meaning of the chiral limit 
in this case is somewhat subtle and this subtlety greatly affects the applicability of our 
analysis. As the current quark masses tend to zero (including m s ) two phenomena happen 
in the nuclear medium. One is the change of properties of particles, which is the subject 
of our analysis. The second phenomenon is the change of the ground state of the medium 
itself. 

In the strict SU(3) chiral limit the ground state of matter has equal amount of up, 
down and strange quarks. In that case the octet of axial vector charges commutes with 
the (baryon number) constraint, and pions, kaons and rj are the usual pseudo-Goldstone 
excitations with E a ~ y/ m q- This is not the situation of relevance at moderately low 
densities. At first sight, one might rather keep the matter non-strange, by imposing the 
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constraint fi u u + u + fidd + d, and watch the formal behavior of the kaon as m s — > 0. In 
this case the a = 4, 7 components of the axial vector charge do not commute with the 
constraint. The situation is analogous to the case of the charged pion in nonsymmetric 
medium. We only have to change the usual I-spin into U-spin or V-spin. Again, only two 
kaon-like excitation (out of 4 in the vacuum) are pseudo-Goldstone modes, with excitation 
energy ~ m q . 

The preceding argument, however, is formal; it depends on using fi u u + u + fidd + d as the 
constraint. This is not the relevant constraint for the neutron matter calculations — at least 
in the interesting applications to neutron stars. In that case of interest the constraints are 
on the net baryon density and the net charge density [p!5H T9| . The reason that these are the 



appropriate constraints is clear — in the neutron star case one is working at long times scales 
for which the system can be expected to be in equilibrium with respect to weak interaction 
processes which can change the flavor quantum numbers. Accordingly, the analysis given 
above does not give any quantitative insight into the problem of kaon condensation. 

In conclusion, we have shown the behavior of pseudo-Goldstone modes in dense isospin- 
asymmetric nuclear matter is very different from the case of the vacuum or isospin symmetric 
matter. Of the three pionic modes which exist in the vacuum only two remain in dense 
isovector matter. Moreover, of these remaining modes, the charged mode is anomalously 
light — having a chiral dimension of one rather than one half. 

This work has been supported by the NSF-Polish Academy of Science grant #INT- 
9313988, NSF PYI grant #PHY-9058487, DOE grant #DE-FG02-40762, and by the Maria 
Sklodowska-Curie grant #PAA/NSF-94-158. 
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